o 
o 
o 

<^ : Interfacial Tension of Electrolyte Solutions 

^ ; Yan Levin 

Institute de Fisica, Universidade Federal do Rio Grande do Sul 

c^ . Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil 

O 
^ February 1, 2008 



o 



'^ . Abstract 

O \ A theory is presented to account for the increase in surface tension 

of water in the presence of electrolyte. Unlike the original "grand- 
canonical" calculation of Onsager and Samaras, which relied on the 
►^ _ Gibbs adsorption isotherm and lead to a result which could only be ex- 

Cn| ■ pressed as an infinite series, our approach is "canonical" and produces 

^^ . an analytic formula for the excess surface tension. For small con- 

Qs^ [ centrations of electrolyte, our result reduces to the Onsager-Samaras 

^^ ' limiting law. 
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1 Introduction 

It has been known for almost a century that aqueous solutions of inorganic 
salts have greater surface tension than pure water [Q . An explanation of this 
effect was advanced by Wagner ||^ in 1924 on the the basis of the theory 
of strong electrolytes, which was introduced only a year earlier by Debye 
and Hiickel . The fundamental insight of Wagner was to realize that the 
presence of ions polarizes the air- water interface, inducing an effective surface 
charge. Since the dielectric constant of water is significantly larger than that 
of air, each ion's image charge equals to it in sign and magnitude. Thus, 
the repulsive interaction with the images reduces the density of electrolyte 
near the interface. Appealing to the Gibbs adsorption isotherm 0, it is 
evident that the depletion of solute near the interface results in an increase 
of surface tension. Furthermore, the experimental observation that for small 
concentrations this increase depends only weakly on the ionic size [|l], §, |^, 
^ ^, suggest the existence of a universal limiting law, similar to the one 
obtained by Debye and Hiickel for bulk properties of electrolyte solutions . 
Indeed, the calculation of Onsager and Samaras {OS) seems to confirm such 
a limiting law 0. A number of approximations adopted by OS in the course 
of calculations, however, obscure the full range of validity of their findings. 
To check the thermodynamic consistency of the 05* results it is, therefore, 
worthwhile to explore other routes to surface tension. In the absence of an 
exact calculation, these will provide a way to asses the self-consistency of the 
formulas obtained. To this end, we propose a new approach for calculating 
the increase in interfacial tension of water due to 1:1 electrolyte. Our method 
differs from that of OS* — who integrated the Gibbs adsorption isotherm — 
in that we identify the excess surface tension directly with the Helmholtz free 
energy necessary to create an interface. The advantage of this approach is 
that it allows to write an analytic formula for the excess surface tension. This 
should be contrasted with the method of OS, who were able to express their 
result only as an infinite sum. Nevertheless, in the limit of large dilution our 
expression reduce to that of OS, suggesting that the limiting law is, indeed, 
exact. The paper is organized as follows: in Sec. H we briefly review the 
thermodynamics of interfaces and the relationship between the "canonical" 
and the "grand-canonical" routes to surface tension; in Sec. |^ we outline 
the OS theory and present our calculations; the results and conclusions are 
summarized in Sec. |[ 



2 Thermodynamics of Interfaces 

Consider an r component mixture confined to volume V. The system is 
in contact with a hypothetical reservoir of solute particles at temperature 
T and chemical potential /ij, with {i = l...r}. If the periodic boundary 
conditions are imposed on V, at equilibrium, the system will have {iVj} 
particles uniformly distributed throughout the volume. 

On the other hand if V forms part of a larger system, its domain must be 
delimited by an interface. We shall idealize this interface as a mathematical 
surface — the Gibbs dividing surface for solvent. The discontinuity produced 
by the interface will affect the interactions between the particles in its vicinity. 
If the interface-particle interaction is repulsive, it will lead to the depletion 
of solute from the interfacial region, forcing it back into reservoir. On the 
other hand, if this interaction is attractive, a concentrated layer of solute 
will build up along the interface, producing a net flow of solute from the 
reservoir. When the equilibrium is reestablished, the distribution of solute 
is no longer uniform, but characterized by a density profile Pi{z), where z is 
a distance from the interface. Clearly if the system has a thermodynamic 
limit, the interface does not influence the bulk distribution of particles and 
Pi(oo) = Pi, where pi = Ni/V. The presence of an interacting interface, 
however, is responsible for a net increase or decrease of solute in the system. 
Thus, we define the amount of solute "adsorbed" as 

n:= / ip.iz)-p.)dz. (1) 



Note that this quantity can be either positive or negative depending on 
whether the solute enters or leaves the system. Now, from general ther- 
modynamic principles, the change in the total Helmholtz free energy of a 
system at fixed volume, temperature, and the amount of solvent is, 

r 

dF = adA + ^Pi dNi, (2) 

where a is the surface tension and A is the area of the interface. We shall 
now divide the solute particles into the "bulk" N^ and the "surface" Nf. 
Eq. H can also be separated into the bulk, 

r 

dF'' = Y,p,dN^, (3) 



and the surface contribution, 

r 

dF' = adA + J2f^idN- . (4) 



i=l 



Since dF'^ is an extensive function of A and {N^}, the Euler's theorem for 
first-order homogeneous functions allows us to integrate Eq. yielding, 



r 



F^ = aA + Y,f^^N!■ (5) 



j=i 



On the other hand, differentiating Eq. ^ and comparing it with Eq. ^ we 
find a Gibbs-Duhem-like equation, 

r 

Ada + ^ A^f d/ii = 0. (6) 

i=l 

For r = 1 this reduces to the Gibbs adsorption isotherm. 
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Since the thermodynamic stability requires dfi/dp > 0, where p is the concen- 
tration of solute, it is evident that a positive adsorption leads to a decrease 
in surface tension, while a negative adsorption increases the surface tension. 
This, then, explains the Wagner's original observation that the repulsion of 
ions from a polarized air-water interface results in a depletion of electrolyte 
and an increase in surface tension. 

Knowledge of Nf allows a calculation of the excess surface tension through 
the integration of Eq. ^. Following Wagner, this was the procedure adopted 
by 05. 

The discussion outlined above relies on the presence of a hypothetical 
reservoir. In the language of statistical mechanics it is intrinsically "grand- 
canonical". A different, "canonical" calculation should also be possible. In 
the thermodynamic limit, the choice of ensemble will not matter, if an exact 
calculation is performed. In practice, however, no exact calculation is possi- 
ble and approximations have to be made. Thus, there is no a priory guarantee 
that the two ensembles will lead to identical results. The canonical approach 



presented bellow, is conceptually simpler than its grand-canonical counter- 
part, since no reservoir is present. Thus, when the interactions between the 
interface and the solute are turned on, the particles do not leave the systems 
and N^ = 0. Therefore, canonically the interfacial tension is equivalent to 
the surface Helmholtz free energy density, a = F'^/A. Now, consider a mix- 
ture confined to a cylinder of length H and a cross-sectional area A. Define 
F as the total Helmholtz free energy of solute and F^"''^ as the free energy 
of solute in the absence of an interface. The change in the surface tension of 
solvent due to addition of solute is then 

a^^^ = lim 4 lim (F - F''"'^) . (8) 

3 The Surface Tension 

We are interested in the surface tension of an interface between an aqueous 
solution of a symmetric 1:1 electrolyte and air. The extension to asymmetric 
electrolytes is, in principle, straight forward. Some extra care, however, has 
to be taken to account for the strong correlations between the cations and 
the anions, which result from an increased ionic charge |1^, 0. In view 



of the experimental observation that for small concentrations of electrolyte 
the excess surface tension depends only weakly on the ionic size ll|, §, §], to 
simplify the calculations we shall treat ions as point-like. The solvent will be 
modeled as a uniform dielectric medium. 

According to statistical mechanics, the concentration of solute a distance 
z from the interface is given by the Boltzmann distribution, 

p,{z) = p,e-^^»(^) , (9) 

where /3 = l/Zc^T, and Wi{z) is the adsorption potential of the specie i and 
is defined as the work required to bring a particle from infinity to distance 
z from the interface. For symmetric electrolyte Wj^{z) = W-{z) = W{z), 
and p+{z) = p-{z) = p{z). Now, consider an ion located at distance z from 
the interface. If the electrolyte is infinitely dilute, the electrostatic potential 
distance ri from the ion can be calculated directly from the Laplace equation. 



where D is the dielectric constant of water, D' is the dielectric constant 
of air, and r2 is the distance from the image charge located opposite to 
the ion at —z. From the second term of Eq. |10|, the charge of the "image 
ion" is qimage = [D — D')q/{D + D'). For an aqueous solution close to 
room temperature D ^ 80 and D' ^ 1, so that qimage ~ q- Therefore, 
in a perturbative theory, D' / D can play a role of a small parameter. Since 
D ^ D' , the zeroth order calculation is already quite accurate, and D' can be 
set to zero. This is the first approximation proposed by Wagner [^ and used 
by OS 0. In the following discussion we shall also adopt this approximation. 
The repulsive force felt by an ion due to the dielectric discontinuity produced 
by an air-water interface is then 

Hz) = T^. (11) 

The amount of work required to bring this ion from infinity to a distance z 
from the interface is 

/z 2 

Hx)dx = ^ , (12) 

where we have added subscript oo to W{z) to emphasize that the calculation 
is done at infinite dilution. Alternatively, the Giintelberg charging process 



12| can be used to calculate the amount of electrostatic work necessary to 



'create" an ion of charge q at distance z from the interface. 



' M 
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Woo{z)= ^qd\ = i-. (13) 



2Dz" ADz 

If the electrolyte is at finite concentration, the electrostatic potential in 
the vicinity of a fixed ion satisfies the Debye-Hiickel equation, 

V^^ = K^^, (14) 

where n'^{z) = 8TTp{z)/DkBT. In order to simplify the calculations Onsager 
and Samaras suggested replacing k{z) in Eq. [l^ by its bulk value /t;(oo) 0. 
This, certainly, seems like a reasonable thing to do in view of the fact that for 
an aqueous solution at room temperature, the boundary layer is very narrow 
and the density profile rapidly approaches its bulk value. Nevertheless, this 
approximation introduces some internal inconsistency into the theory which 
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can, in principle, manifest itself when different thermodynamic routes are 
taken to calculate the excess surface tension. To zeroth order in D' / D and 
with k{z) — * fi:(oo), Eq. O can be integrated pi, |13| yielding 



V^(ri,r,) = -^^ + -^^, (15) 

where, once again, ri is the distance from the ion, and r2 is the distance from 
the image charge located at —z. The adsorption potential is obtained through 
the Giintelberg charging process or by direct force integration producing 

W(.) = ij^. (16) 

The density profile for cations and anions is found by substituting this ex- 
pression into Eq. ^ 

P±(^) =pexp ( ^— j . (17) 

We have defined h as half the Bjerrum length, h = Ab/2 = q^/2DkBT. For 
water at room temperature b ~ 3.6 A — comparable to the size of a hydrated 
ion — and the density profile rapidly reaches its bulk value. The amount 
of solute adsorbed can now be calculated by inserting Eq. |1^ into Eq. |l|. 
With skillful complex variable analysis, OS were able to integrate the Gibbs 
adsorption isotherm, obtaining an expression for the excess surface tension 
as an infinite series in /tfe 0. As was discussed in Sec. ||, this procedure 
does not conserve the number of particles in the system and is, intrinsically, 
grand-canonical. We now present an alternative, canonical, calculation of 
the excess surface tension. 

Lets suppose that the electrolyte is confined to a cylinder of height H and 
cross-sectional area A. If the interface-ion interactions are neglected (periodic 
boundary conditions), the electrolyte will be uniformly distributed over the 
volume of the cylinder. In the thermodynamic limit, the electrostatic free 
energy can be easily calculated from the Debye-Hiickel theory yielding, 

2 

where A''^ = A^+ + A^_ is the total number of solute particles. On the other 
hand, the presence of an interface produces a concentration gradient charac- 



terized by the normalized distribution, 

Now, suppose we fix an ion some distance z from tlie interface, far from the 
radial boundary of the cyhnder. The electrostatic potential in the vicinity 
of this ion is approximately given by Eq. ^ — where in view of the thermo- 
dynamic limit of Eq. p, we can neglect the finite size corrections. Evidently 
this potential is produced by the ion itself, as well as, by the interface and 
the ionic atmosphere. The potential that the ion feels due to polarization of 
the interface and the ionic atmosphere is, 

V.„W=.im(v-(n.--.)-^)=^ + ^. (20) 

The electrostatic energy is, 

^ = ^\ {p^{z)+p_{z))M^)dz. (21) 

This can be subdivided into the bulk and surface contributions, correspond- 
ing to the first and the second term of Eq. |20|, respectively. The bulk term 
is easily integrated yielding, 

E^ = ^iV,. (22) 

2D ^ ' 

The surface contribution is found to be. 



f" dz 

, / exp(-2fi:z-/?iy(z)) — 



E' = ^—^^ rr . (23) 

AD f" ^ ' 

/ e-x.Y>{-(3W{z))dz 
Jo 

Now, in the limit H ^ oo relevant for the calculation of surface tension, 

lim -- / (e-^^(^) -l)dz = 0, (24) 



and the Eq. ^ simplifies to, 

Unfortunately, no method is available for analytic evaluation of this integral. 
We note, however, that if the adsorption potential is replaced by its value at 
infinite dilution, W{z) — > Woo{z), the integral can be done explicitly yielding, 

E' = ^Ko (2v^) , (26) 

where K^{x) is the modified Bessel function of order u. Replacement of 
the adsorption potential by its value at infinite dilution should be a good 
approximation, since the large distances — for which the discrepancy between 
W{z) and Woo{z) becomes significant — are not important because of the 
exponential drop in the electrostatic potential away from the interface. To 
confirm this, we have numerically evaluated the integral 

I(Kb)= re-2--^^(^)— (27) 



and compared it with the exact analytic expression obtained when W{z) — »• 
Woo{z), Fig. 1. As was hoped the agreement is, indeed, quite good extending 
all the way to y = nb ^ .45. In water at room temperature, this corresponds 
to concentration of 0.15M, which is above the maximum for which the lim- 
iting laws of this paper can be realistically expected to apply. 

The electrostatic free energy can now be obtained through the Debye 
charging process, in which all the particles are simultaneously charged from 
zero to their full charge]^, [T0[| , 



\E{\q)^, (28) 

note that both b and n are dependent on q, so that b{\q) = \^b{q) and 
n{\q) = XK{q). The integral can, once again, be done explicitly yielding 
the electrostatic free energy of an electrolyte solution in the presence of an 
interface, 

q^K^^ - 4^ [2M2^) ,F,{1- 2/3, 5/3; y) + 



3D 2D 

3^K,i2^ iF2(l; 5/3, 5/3; y)] . (29) 
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Figure 1: The solid curve is the numerically evaluated integral I{y), as a 
function oi y = nb, Eq. ^; the dashed curve is the analytic expression 
obtained when W{z) — > Woo{z), see Eq. |20 



Here, pFg is the generalized hypergeometric function [jT4| and y = bn. Sub- 
stituting Eq. ^ into Eq. ||, the increase in surface tension of water due to 
1:1 electrolyte is, 

cr'^^ = < [2Ko{2^) iF2(l; 2/3, 5/3; y) + 3^K,i2^) 1^2(1; 5/3, 5/3; y)] ,(30) 

where a^^ = q'^p/2D. For very low concentrations, this expression reduces to 
the limiting law found by OS, 



where •je = 0.57721566490153... is the Euler's constant. 



(31) 



4 Results and Discussion 

We have presented a "canonical" calculation of the excess surface tension 
in an electrolyte solution. Unlike the earlier "grand-canonical" method of 
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Onsager and Samaras, our approach leads to an analytic expression for the 
excess surface tension expressed in terms of Bessel and hypergeometric func- 
tions. It is gratifying, however, that in spite of all the approximations, the 
two ensembles produce the identical limiting law. This thermodynamic self- 
consistency suggests that the OS limiting law is, indeed, exact to zeroth 
order in D' / D. In Fig. 2 we plot, as a function of concentration, the a^^ 
Eq. |30|; the sum of the first twenty terms of the infinite series for the excess 
surface tension obtained by OS 0; and the OS limiting law Eq. |31|. The 
surface tension a is measured in mN ■ m~^ and the concentration of salt c in 
moles/liter (M) so that, 

69.4692c. ^^ _i. 

For water at room temperature D ^ 78.54. We note that both ensembles 
agree fairly well over the full range of concentrations, with the "canonical" 
calculation predicting a somewhat larger excess surface tension. In Fig. 1 it 
was shown that the substitution W{z) -^ Woo{z) leads to a good approxi- 
mation for the electrostatic energy. To confirm that this also extends to the 
free energy, we have performed the Debye charging process, Eq. ^, numeri- 
cally using the explicit form of W{z), Eq. |T6|. The result is plotted as solid 
points in Fig. 2. Indeed, the numerically calculated surface tension agrees 
quite well with the analytic result, Eq. ^, obtained with the substitution 
W{z) ^ W^{z). 

It has been noted that the 05* theory gives a fairly good quantitative 
description for concentrations up to O.IM, above which it consistently un- 
derestimates the increase in the interfacial tension [^ |^, ^, Q. The canonical 
calculation presented above extends the range of agreement between theory 
and experiment. It is, however, unrealistic to demand that the theory pre- 
sented above should apply to concentrated solutions, for which even the bulk 
thermodynamic properties loose their universality. Thus, for concentrations 
above 0.2M, the molecular nature of the solvent as well as the lyotropic 
properties of solute will become important. In fact, it has been observed 
experimentally that for concentrated solutions the excess surface tension in- 
creases linearlywith the concentration of electrolyte. For an aqueous solution 
of NaCl at T = 25°C it is found that a^^ ~ 1.6c [jT5[, which begins to show 
strong deviation from Eq. ^ for c > 0.2M. 
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Figure 2: The excess surface tension of an aqueous solution of a 1:1 electrolyte 
at room temperature. The solid curve is the analytic expression given by 
Eq. 13^; the long dashed curve is the sum of the first 20 terms of the OS result 
0; the short dashed curve is the Onsager-Samaras limiting law, Eq. |^; the 
solid points are the result of numerical integration of Eq. ^ with full W{z). 
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6 Figure Captions 

Fig. 1: The solid curve is the numerically evaluated integral I{y), as a 
function oi y = nb, Eq. ^; the dashed curve is the analytic expression 
obtained when W{z) -^ Woo{z), see Eq. |2^ 

Fig. 2: The excess surface tension of an aqueous solution of a 1:1 elec- 
trolyte at room temperature. The solid curve is the analytic expression given 
by Eq. ^ the long dashed curve is the sum of the first 20 terms of the OS re- 
sult 0; the short dashed curve is the Onsager- Samaras limiting law, Eq. |31|; 
the solid points are the result of numerical integration of Eq. |2^ with full 
W{z). 
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